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Abstract. In this paper we prove that weak solution constructed by artificial compress- 
ibility method are suitable in the sense of Scheffer, [T7], [18]. Using Hilbertian setting and 
Fourier transform with respect to the time we obtain nontrivial estimates of the pressure 
and the time derivate which allow us to pass into the limit. 



1. Introduction 

In this paper we investigate if the weak solutions of the Navier-Stokes equations are suitable 
in the sense of Scheffer [T7], [18]. The incompressible Navier-Stokes equations in three spatial 
dimensions with unit viscosity and zero external force are given by the following system 

idtu- Au + {u-\7)u + \7p = ^^^^ 
] div M = 0, 

where {x,t) G M'^ x [0,T], u £ M'^ denotes the velocity vector field and p G M the pressure of 
the fluid. Let us recall the notion of Leray weak solution. 

Definition 1.1. We say that u G L°°{{0,T); L'^{R^)) n L'^{{0,T); H'^{R^)) is a Leray weak 
solution of the Navier-Stokes equations if it satisfies (|1.1|) in the sense of distribution for all 
£ (M^ X M) with div ip = and moreover the following energy inequality holds for every 
t£[0,T] 

dx\u{x,t)\'^ + 2 ds dx\Vu{x,s)\'^ < dx|ii(x,0)|^ (1.2) 

Jo Jk3 

In the mathematical literature there exists several proofs of the global existence of Leray 
weak solution for divergence- free initial data in L^(R^), see for example the book of P. L. Lions 
|13j or the monograph of Temam [26j. Several problems about the Leray weak solutions are 
still open, for example it is not know whether or not the solutions are unique or develop 
singularities in a finite time also for smooth initial data. In fact the regularity requirements 
for proving uniqueness and global in time regularity are not available at the present time for 
Leray weak solutions. Scheffer in pT], [18] introduced the notion of suitable weak solutions 
that we recall here. 

Definition 1.2. Let (n,p), u G L'^{{0,T); H'^{R^))nL'^{{0,T); L'^{R^)), p £ V' {{0,T); L^{R^)), 
be a weak solution to the Navier-Stokes equation (jl.ip . The pair {u,p) is said a suitable weak 
solutions if the following local energy balance 

dtl^luA +V - ((^\u\^+p\u] -A(^\uA+\Vu\^-f-u<0 (1.3) 
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holds in the distributional sense. 

In mathematical literature (|1.3p is called generalized energy inequality. 
It is important to remark that weak solutions constructed by Leray are suitable. At the mo- 
ment the best regularity result for the weak solution of the Navier-Stokes equations is a partial 
regularity result, i.e. the so called Caffarelli-Kohn-Nirenberg theorem, |2j. This theorem as- 
serts that the one-dimensional parabolic HausdorfF measure of the singular set is zero and 
holds only for suitable weak solutions. If the class of suitable weak solutions is a proper class 
of Leray weak solution is an open problem since Scheffer's works. The method usually used 
to construct suitable weak solutions are regularization of the non- linear term , [2], adding 
hyper viscosity, [1]. Recently Guermond in [8] proved that some type of Faedo-Galerkin 
approximation lead to a suitable weak solution. He used classical Lions method of frac- 
tional derivate in order to obtain maximal regularity in negative Sobolev space. In particular 
Guermond proved that Faedo-Galerkin weak solutions of the three dimensional Navier-Stokes 
equations with Dirichlet boundary condition are suitable provided they are constructed using 
finite-dimensional approximation spaces having a discrete commutator property. 

In this paper we prove that Leray weak solutions costructed by the artificial compress- 
ibility method are suitable in the sense of the Definition 1.2. The artificial compressibility 
approximation was introduced by Chorin [31 U], Oskolkov [T5] and Temam |24[ I25| . in or- 
der to deal with the difficulty induced by the incompressibility constraints in the numerical 
approximation. The approximation system reads as follows 

{dtu^ + V/ = Au^ - (n^ • V) - ;^(divii^)ti'' 
2 (1-4) 
edtp'^ + divn"^ = 0, 

where {x,t) G x [0,T], = u''{x,t) G and = p^{x,t) eR, f = G M^. 

Temam showed the convergence of this approximation on bounded domain. Recently in [6] and 
[7] the result was extended in the case of the whole space and exterior domain, respectively. In 
[6] the convergence towards Leray weak solutions of the Navier-Stokes equations is achieved by 
using the dispersive structure of the system. Here we will consider the system ()1.4p endowed 
with the following initial conditions 

u'ix,0)=ulix) p'ix,0)=plix), (1.5) 

such that 

ul uo in L^{R^) as e ^ 0, (1.6) 

^/epl^Om L^{R) as e ^ 0. (1.7) 

As we will see later on in Section 3, in order to get some a priori estimate on and we 
need to assume that 

G H\R^). (1.8) 
We will be able to prove the following theorem 

Theorem 1.3. Let {u^,p^) he a weak solution of the system ()1.4p with initial data (jl.Sp 

and such that ()1.6p . ()1.7p and ()1.8p are satisfied. Then {u'^ ,p'^) converges to a suitable weak 
solution of the Navier-Stokes system as e goes to zero. 

In order to prove this theorem we have to estimate carefully the pressure term p^ and yfep^ ■ 
Since we haven't the incompressibility constraint we cannot use classical method based on the 
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elliptic equation associated to the pressure. Also the dispersive approach, as in [6j, doesn't 
give usefull estimates. We will use the method of Lions of the fractional derivates for getting 
the necessary estimates. This paper is organized as follows. In Section 2 we recall some basic 
facts about Navier-Stokes equations and the artificial compressibility method and fix some 
notations. In Section 3 we obtain the a priori estimates which allow us to pass into the limit. 
In Section 4 we give the proof of the main result. We want to point out that the estimates of 
the Lemma 3.4 are the core of this paper. 

2. Notations and preliminary results 

For convenience of the reader we establish some notations and recall some basic fact about 
the Navier-Stokes equations. 

2.1. Notations. We will denote by V{W^ x M+) the space of test function C^{W^ x M+), by 
V'iR"^ X M+) the space of Schwartz distributions and (•, •) the duality bracket between D' and 
V and by MfX' the space C°([0, T]; X)'. The inner product in L^{W^) wih be denoted by 
parentheses, namely {u,v) := f^^u{x)v(x)dx. Moreover W'^'P{M.'^) = (/ — A)~2 LP(M°') and 
H''{R'^) = VF*^'2(M<^) denote the non homogeneous Sobolev spaces for any 1 < p < oo and 

fc G R. W^'PiR"^) = (-A)tLP(R'^) and H''{R'^) = W''''^{R'^) denote the homogeneous Sobolev 
spaces. 

Let -ff be a Hilbert space, 6 > 1, L^(R;H) is the Lebesgue space of function with values in 
H, namely 

L\R;H) := ■.R3t^i;{t) e H; j ||V'(t)||f^dt < ooj . 

For any •i/^ G L-'^(M; H) and A; G M we define the Fourier transform of il) with respect to time 
as follows 

Jr 

The previous definition can be extended to the space of tempered distribution S (R; H) in 
the usual way. We define H'^(R;H) the space of tempered distributions v £ S {R; H) such 
that 

/ il + \k\)^^\\v\\jjdk<+oo. 
Jr 

The space H'^{{0,T); H) is defined by those distributions that can be extended to S (R;H) 
and whose extension is in H^{{0,T);H). The norm in H"' ({0,T); H) is the quotient norm, 
namely, 

\\v\\H-'(iO,T);H) = inf \M\H'y{R;H)- 

v=u 

q.o. su (0,T) 

The notations LP{L'^), LP{W^''^), Lf{H^) and H'''{H^) will abbreviate respectively the spaces 
LP{[0,T];Li{R'^)), LP{[0,T];W''''i{R'^)), LP{[0,T]; {R'^)) and H'' {[0,T]; (R'^)). Moreover, 
we shall denote by Q and P respectively the Leray's projectors Q on the space of gradients 
vector fields and P on the space of divergence - free vector fields. Namely 

Q = VA-Miv P = I-Q. (2.1) 

Let us remark that Q and P can be expressed in terms of Riesz multipliers, therefore they 
are bounded linear operators on every W^'P {1 < p < oo) space (see p3]). 
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2.2. Compactness lemma. We recall some Aubin-Lions like compactness results. We will 
use these lemmas in the proof of Theorem 1.1. 

Lemma 2.1. Let Hq C H C Hi be three Hilbert space with dense and continuos embedding. 
Assume the the embedding Hq G H is compact and let ^ > be a positive real number. Then, 
the injection L^{{0,T); Hq) n H^{{0,T); Hi) L^{{0,T);H) is compact. 

Proof. See Lions [9j p. 61 Theorem 5.2 □ 

Lemma 2.2. Let X C Y be two Hilbert space with compact embedding and let t > ^. 
The injection H'^{{0,T); X) — t- C^{[0,T]; X) is continuos and the injection H'^ {{0,T); X) — )■ 
C^{[0,T];Y) is compact. 

Proof. See |8], Appendix A.l. □ 

Lemma 2.3. Let Hq C Hi be two Hilbert spaces with compact embedding. Let 7 > and 
^ > u, then the injection H"'{{0,T); Hq) C H'^{{0,T); Hi) is compact. 

Proof See [8j Appendix A.2 □ 

2.3. Artificial Compressibility Approximation. In this section we recall some previous 
result about the approximating system that we rewrite for convenience of the reader. 

5tn" + V/ = fxAu"" - {u" • V) - ^(div 2) 
edtp^ + div = 0, 
As said in the Introduction we consider two initial condition, namely 

u'{x,0)=u'q{x) p'{x,0)=p'q{x), (2.3) 

such that 

uf, uq in l2(M^) as e (2.4) 

VePo ^ in L'^(R) as e ^ 0. (2.5) 

We have an existence theorem for the system ()2.2p with the above initial data. 

Theorem 2.4. Let {uq,Pq) satisfy the conditions (|2.4p and (j2.5|) . e > 0. Then the system 
(12. 2p has a weak solution (n^jp"^) with the following properties 

• e L~((0,T);L2(r3)) nL2((0,r);/i'HR3)) 

• VeP^ G L°°((0,r);L2(M3)) 

for all T >0. 

The proof of this theorem can be omitted since it is a direct consequence of all the a priori 
bounds in fB] and it follows by using standard finite dimensional Galerkin type approximations. 
In [B] the following result was showed. 

Theorem 2.5. Let {u^,p^) be a sequence of weak solutions in M'^ of the system ()2.2p . assume 
that the initial data satisfy (|2.4p and ()2.5p . Then 

• There exist u G L°°{{0,T); L^{R^)) n L^{{0,T); H^{R^)) such that 

u" weakly in L'^{{0,T); H^{R^)). (2.6) 
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• For any p £ [4, 6) the gradient component Qu^ of the vector field u'^ satisfies 

Qu" strongly in L^{{0,T); LP{R^)). (2.7) 

• The divergence-free component Pu^ of the vector field satisfies 

Pu' ^ Pu = u strongly in L'^{{0,T); LI^{R^)). (2.8) 

• The sequence {p"^} converges in the sense of the distributions to 

p = A-^tr {{Duf}) (2.9) 
Moreover Pu = u is a Leray weak solution to the incompressible Navier-Stokes equation 

P {dtu -Au + {u- V)m) = in L>'([0, T] x R^) (2.10) 
and the following inequality holds for every t £ [0,T], 

dx\u{x,t)\^ + 2 [ ds [ dx\Vu{x,s)\^ < [ dx\u{x,0)\^. (2.11) 
Jo Jr3 Jm.3 

It is worth to mention here that in order to prove Theorem 2.4 and Theorem 2.5 only the 

conditions ()2.4p and ()2.5p are needed. We don't need to assume Uq S H^, namely (jl.Sp . 

3. A Priori Estimates 

In this section we deal with the a priori bounds for our approximating system. 
By taking into account the Corollary 4.2 and Theorem 4.3 of [6j we recall the following 
estimates. 

Lemma 3.1. Let (u'^p^) be a weak solution of the system ()2.2p with initial data (|2.3p . Assume 
that the conditions (|2.4p and (12. 5p hold. Then there exists c > 0, indipendent on e, such that 
the follwing estimates hold 





is bounded in L°°([0, T]; L2(I 






is relatively compact in /7^^([0, T] x 


R^), 




is bounded in L^([0, T] x 






is bounded in L°°{[0,T]; L'^{R^)) n L'^{[0,T]; L^(I 


^')), 




is bounded in ( [0, T] ; (M^) ) n ( [0, T] ; L^/^ 


^')), 


{divu'^X 


is bounded in ( [0, T] ; (M^) ) n ( [0, T] ; L^/^ (E 






is bounded in L'^{[0,T];W-'^'^{ 




e^l^dtp' 


is bounded in L^([0, T]; I^"^'^( 





Unfortunately these estimates are not sufficient to prove the generalized energy inequality 
(|1.3p . The aim of this section is to estimate more carefully the nonlinear term and to get 
estimates on Hilbert space for the pressure p^ and the term \fep^ . 
Let be p^ g, f and s real numbers such that the following relations holds: 

2 3 , , , 3, s 1 1 11 , , 

- + - = 4, p. [1.21. .S^H.^I. 3---2. ''-^-2- 



These relations follow from the Sobolev embedding theorems, in particular if p, q satisfy ()3.1 
then the following embedding holds 

LP((0,r);L'?(M3)) c i7-'^((0,r);F-"(M=^)) (3.2 
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The first a priori estimate regards the nonlinear terms 

Lemma 3.2. Let he a weak solutions of the system ()2.2p . then there exists a constat c > 0, 
indipendent on e, such that 

\\{U' ■V)U' + ]^U' diYu'\\H-ri^H-s-^<C. (3.3) 

The proof of ()3.3p can be omitted since it follows by standard interpolation argument and 
the embedding ()3.2p . 



3.1. Estimates of the velocity vector field. In this section we get the estimate on and 
(9tn^ on negative Sobolev spaces. So we need to assume 

Since we are going to use Fourier transform with respect to time we need to extend all the 
function from [0, T] to M. We define by u'^ the following extension of 

(t+lK on [-1,0] 

(o,r + i) 

on[r + l,oo]. 

Let if G C°°(M) be such that supp{ip) C {—1,T + 1) and if = 1 on [0,T], we denote with 
abuse of notation 



Next we define the following function 



fs _ i (l + t)y?'u§ + v9n§-(l + t)v7((n^-V)u^ + iM^divw'^) t € {-1,0) 
\ -93((n^-V)n^ + iu^divn^) + (^'u^ t ^ (-1,0) 

It follows that e are well definied on (— oo,+oo). By using (jl.Sp it follows that there 
exist c > 0, indipendent on e, such that 

<c (3.4) 

Lemma 3.3. Let {u^ ,p^) be a solution of the system ()2.2p then, there exists c > 0, independent 
on e, such that 

• for all s G [^) f ) ^^'^ r > f, one has 

||<9t'u^||j^^,(^_,) < c, (3.5) 
\\^'^^\\h-{h-=) ^ (3.6) 

• for all a G \_\^\) o-nd r < f = |(1 + a), one has 

Proof. We start by proving ()3.7p . By taking the Fourier transform with respect to time of 
the system ()2.2p we obtain 



jlinku^ - Aii^ + Vp^ = 

< ~ ~ (3.8) 

I elinkp'^ + div u'^ = 



ARTIFICIAL COMPRESSIBILITY AND SUITABLE WEAK SOLUTIONS 



Let a > 0, we multiply the first equation of ()3.8p by the complex conjugate of — A~"m^ 
and the second by the complex conjugate of — A~°p^. By summing up and by taking the 
imaginary part we obtain the following inequality 

|A;||(n^A-"n-)| < |(^, A-%~-)|2. (3.9) 

Now we choose a such that a < s < 1 + 2a. Since s G [|, |) we have that a G [|, By 
interpolation we get 

lln-ll^. <||n-|i;^Jk~1l^-T,<., (3.10) 
where 7 = Inserting ([3?T0|) in (f3r9|) we have 

|fc|||n-||?r-"^<c||/,||^,-.||n~^||^-,^. (3.11) 

We set v = TT"^ , then we have 

\k\^-''\\u-\\l-. < c(l + |A^|)-1lMl^J|n-||^. (3.12) 
By integrating (|3.12|) with respect to the time and by using Holder inequality we get 

l^dk\k\^\\u^\\l.^ < ciirii;p,(^_.)iin-ii;t^). (3.13) 

If we set r := i±f (1 - f) then ([SZD follows. 

Now we are going to prove ()3.5p . We multiply the equations of the system ()3.8p respectively by 
the complex conjugate of — A~''m^ and — A~'*p^. By summing up and by taking the imaginary 
part we obtain 

|A;|||n-|li-. <c||f ||^-.||A-^n-||^.. 

So we have 



\\dM\H'S<c\\f\\H-s. (3.14) 



Thus ()3.5p follows by integrating with respect to time (|3.14p . It remains only to prove (|3.6p . 
Testing the first and the second equation of ()3.8p with —A^^'^u^ and with A^~'^p^ respectively, 
we have 

\\Au'\\h-. < c||/^||j^-.||Au^||h— (3.15) 
Thus by integrating in time (j3.15p we obtain (|3.6p . □ 

3.2. Estimates of the pressure. 

Lemma 3.4. Let {u^ ,p^) be a weak solution of (|2.2p . there exist a constant c > 0, independent 
on e, such that for s E |), r > f = | — |, 5 G (0, ^) and /3 G (O, ^—^^) we have 

MWH-rim-") ^ c (3.16) 

and 

■ 1 + < c (3.17) 
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Proof. The estimate ()3.16p can be proved by observing that from the first equation of the 
system ()2.2p we have 

V/ = f + Au" - dtu"" (3.18) 

and by using (j3.5p , (j3.6p and ()3.4p . Now we are going to prove (j3.17p . We have the fohowing 
equation 

2mkep^ + div -u^ = (3.19) 
By multiplying (|3.19p by A^a+^p"^ we obtain 

\mV~ep^\t.l,^, < c|(A-HV, V/)|. (3.20) 



From (j3.20p we have 



\\P' 



1 S 



\k\'+h\V~ep%.^, < (1 + \k\)h+^\m ^ fil^, (3.21) 

^ ^ " (l + |A:|)2 + 2 

where ^ G (O, ^—^)- By integrating with respect to time ()3.2ip and by using the Holder 
inequality we have 

The right hand-side of ()3.22p is bounded by using (|3:7P and (fXT6|l . provided (5 G (0, ^) and 

/3e(o,i^). □ 

4. Convergence to a suitable weak solution 

In this section we give the proof of the theorem 1.1. 
Let us multiply the first equation of (j2.2p by u'^cp with (p G Co°(M^ x E), (/> > 0, then we have, 

r 1 

/ dt{dtu^,u''(l)) - (An^^i^0) + (n^ • V^i^^i^0) + -(n= div n^ u^c/)) + (Vp^u^(^) = 0. (4.1) 
Jo 2 

By integrating by parts we obtain 

'Z"* I 1 2 I 1 2 

/ dt(|Vn^|2,</.) = [ d^^(</.t + A,/.) + (^/^^,V0) + (nV,V0) + (/divn^0). (4.2) 
We estimate each term of ()4.2p separately. 

By weak lower semicontinuity and the fact that ^ u weakly in L^{H^) we have that 

[ dtdV'up,^) < liminf /" dt(|Vu^p,(/.) (4.3) 
Jo Jo 

Since ^ u strongly in L'^{Lf^^), we get 

^^((/.t + A0)^ / dt^((/.t + A(/.) ase^O. (4.4) 



^ Jo 

2/ r3 



Next, by interpolation we have that -u^ — t- n strongly in L (L ) and that is bounded in 
L^{L^), so it follows 

/ dt('u^-^ — / dt('u— ,(/>) ase^O. (4.5) 
JO 2 Jq 2 
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In order to estimate the last two terms in ()4.2p we have to use carefully the estimates of the 
Lemma 3.3 and Lemma 3.4. We start by estimating 

Jo 

Let r/ > 0, we set r = |+?7, s = ^. This choice implies that is uniformly bounded. 

So ^ p weakly in R-'' {{0,T), H^^^(M.^)). Now let rj' £ [0, ^] and set a = ^ - r/' < i, 
T = |(1 + a — e ). This choice implies that Hii^Hz/r^ji^-ti) is uniformly bounded and as a 
consequence ^ u weakly in H'^{{0, T); H~"). By using Lemma 2.3 we obtain that u'^ ^ u 
strongly in H^'{{0,T); Hj^^{M.^)) provided r > r and s > a. Since t] G (0, ^) we have s > a, 
and if we choose rj = ^(3 — 20r/') we have t > r. Now it follows easily that 



T 



T 



dt{pu,V(f)) 



as e 



0. 



(4.6) 



By using the second equation of ()2.2p we have that 



dt{p'^ div u"^ , ( 



<V~e 


L 







rT 



dt{yfep'^ ,p^ 



<^e\\^ep' 



So by choosing 5 and j3 as in the Lemma 3.4 we have that 

l-T 

dt{p'^ div M^, 0) — )- as e 



0. 
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